Computational Study of Turbulent-Laminar Patterns in Couette Flow 
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Turbulent-laminar patterns near transition are simulated in plane Couette flow using an exten- 
sion of the minimal flow unit methodology. Computational domains are of minimal size in two 
directions but large in the third. The long direction can be tilted at any prescribed angle to the 
streamwise direction. Three types of patterned states are found and studied: periodic, localized, 
and intermittent. These correspond closely to observations in large aspect ratio experiments. 
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Plane Couette flow - the flow between two infinite par- 
allel plates moving in opposite directions - undergoes a 
subcritical (discontinuous) transition from laminar flow 
to turbulence as the Reynolds number is increased. Due 
to its simplicity, this flow has long served as one of the 
canonical examples for understanding shear turbulence 
and the subcritical transition process typical of channel 

and pipe flows [IQ, 111! IIS HUH E3,[ilE3 ° nl y 

recently was it discovered in very lar ge a spect ratio ex- 
periments by Prigent et al. 0, IdH llflj that this flow 
also exhibits remarkable pattern formation near transi- 
tion. Figure ^ shows such a pattern, not from experi- 
ment, but from numerical computations reported here. 
An essentially steady, spatially periodic pattern of dis- 
tinct regions of turbulent and laminar flow emerges spon- 
taneously from uniform turbulence as the Reynolds num- 
ber is decreased. It now appears that turbulent-laminar 
patterns are inevitable intermediate states on the route 
from turbulent to laminar flow in large aspect ratio plane 
Couette flow. 

Related patterns have a long history in fluid dynamics. 
In Taylor-Couette flow between counter-rotating cylin- 
ders, Coles E3 first discovered a state known as spi- 
ral turbulence with coexisting turbulent and laminar re- 
gions. This state was famously commented on by Feyn- 
man [17j and has attracted attention as an example of a 
coherent structure comprising both turbulence and long- 
range order [l8|, UJJ |2CJ, |2JJ . Until recently all experimen- 
tal studies of this state showed only one turbulent and 
one laminar patch. Prigent et al. |l5j found that 

in a very large-aspect-ratio Taylor-Couette system, the 
turbulent and laminar regions form a periodic pattern, 
of which the original observations of Coles comprised 
only one wavelength. Cros and Le Gal |22j discovered 
large-scale turbulent spirals as well, in experiments on 
the shear flow between a stationary and a rotating disk. 
The Reynolds-number thresholds, wavelengths, and an- 
gles are very similar for all of these turbulent patterned 
flows. Moreover^ Prigent et al. suggest that the turbu- 
lent spots [1 II II II S El El El E3 long known to exist 
near transition are essentially a manifestation of the same 
mechanism. 
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FIG. 1: Turbulent-laminar pattern at Reynolds number 350. 
The computational domain (outlined in white, aligned along 
x', z') is repeated periodically to tile an extended region. The 
kinetic energy is visualized in a plane midway between and 
parallel to the plates moving in the streamwise (x) direction. 
Uniform gray or blue corresponds to laminar flow. The sides 
of the image are 60 times the plate separation L y = 2; the 
pattern wavelength is 20 L y . Streamwise streaks (on the scale 
of L y ) are visible at the edges of the turbulent regions. 



We report the first direct numerical simulation of 
turbulent-laminar patterns. Our simulations are de- 
signed to reduce computational expense, to establish 
minimal physical conditions necessary to produce these 
large-scale patterns, and to impose and thereby investi- 
gate the pattern wavelength and orientation. 

Our study extends minimal-flow-unit (MFU) simula- 
tions of turbulent channel flows El & E1 and we begin by 
recalling these. The plates located at y = ±1 move at 
unit velocities ±x. The Reynolds number is Re = l/v, 
where v is the kinematic viscosity of the fluid. The sim- 
ple Couette solution uc = yx is linearly stable for all 
values of Re. However, above a critical Re near 325 [||, 
transition to turbulence occurs for sufficiently large per- 
turbations. The turbulence is characterized by the cycli- 
cal generation and breakdown of streaks by streamwise- 
oriented vortices with a natural spanwise pair spacing 
of about 4 0, H, IS HI El. In the MFU approach, a 



2 



X (uLreamwist;) 




FIG. 2: Simulation domains. The wall-normal direction y is 
not seen; L y = 2. The gray or colored bars represent stream- 
wise vortex pairs with a spanwise spacing of 4. (The vortices 
are schematic; these are dynamic features of the actual flow.) 
(a) MFU domain of size 10 x 4. (b) Central portion of a 
domain [on the same scale as (a)] titled to the streamwise 
direction, a, a 1 and /3, j3' are pairs of points identified under 
periodic boundary conditions in x' . (c) Full tilted domain 
with L x i = 10, L z , = 120, 9 = 24°. On this scale the MFU 
domain, shown for comparison, is small. 

periodic domain of minimal lateral dimensions is sought 
which can sustain this basic cycle. For plane Couette 
flow near transition, the currently accepted 9] size is ap- 
proximately L x x L z = 10 x 4 [Fig. [2(a)]. 

We extend the MFU computations in two ways. First 
we tilt the domain at angle 9 to the streamwise direction 
[Fig- 13b)]. We designate by x' and z' the periodic direc- 
tions of the tilted domain. To respect the spanwise streak 
spacing while imposing periodic boundary conditions in 
x', the domain satisfies L x /sm9 ~ 4, for 9 > 0. (For 
9 = 0, we only require L x i ^ 10.) Secondly, we greatly 
extend one of the dimensions past the MFU requirement 
[Fig. 12(c)] . In practice we use L z / between 30 and 220, 
usually 120. We can thus simulate large length scales 
oblique to the streamwise direction. 

The incompressible Navier- Stokes equations are sim- 
ulated using a spectral-element (x'-y) - Fourier {z') 
code [2^]. The boundary conditions are no-slip at the 
moving plates and periodic in the x' and z directions. 
The spatial resolution for the L x i xL y x L z > — 10 x 2 x 120 
domain is N x x N y x N z — 61 x 31 x 1024, consistent with 
that previously used at these values of Re d, S @ ■ 
suits have also been verfied at higher resolutions. 

We make two comments distinguishing our approach. 
Experimentalists [l3l flil Hif varied Re and reported the 
properties of the resulting patterns: in particular they 
measured angles and wavelengths varying from 9 = 25° 
and \ 7 j = 46 at Re = 394 to 9 = 37° and \ 7j = 60 at 
Re = 340. (They extrapolated the domain of existence 
to be 325 < Re < 415.) In contrast, we fix the pattern 
angle and wavelength: in this way, we can determine the 
boundaries in parameter space within which each pattern 
can exist. Second, all the turbulent states we report are 
bistable with simple Couette flow. A major goal 0, 0, 
[ToL Hlj. not addressed here, has been the determination 



of lifetimes and transition probabilities of turbulent flow 
as a function of amplitude and Re. 

We begin with simulations exploring the dependence 
of patterns on Re. To allow the system sufficient free- 
dom to select different states, we set L z > — 120, two to 
three times the experimentally observed wavelength. We 
fix 9 = 24°, near its observed value at pattern onset. 
Figures EJa) and (b) show two long scries of simulations 
spanning the range 290 < Re < 500. Space-time di- 
agrams are shown for decreasing and increasing Re in 
discrete steps over time. In each case, kinetic energy 
fluctuations are on the right and principle peaks in asso- 
ciated spatial Fourier transforms are on the left. 

More specifically, we compute E = |u — uc| 2 /2 at 32 
points equally spaced in z' along a line [x 1 = y = 0) in 
the mid-channel. We compute E rms , the rms of E in time 
windows of size T = 250. This gives a measure of the 
flow's turbulent intensity on a space-time grid. (Other 
measures such as the rms of individual velocity compo- 
nents gives similar results.) Time windows in Fig. CJa) 
show E from which E ms is computed at two points on 
the space-time grid. For the spectra on the left, we com- 
pute the instantaneous spatial Fourier transform E m of 
E for the same 32 points in the mid-channel. We take 
the modulus (to remove phase information) and average 
over windows of length T — 500 to obtain (\E m \). 

In Fig. 12(a) a turbulent flow is initialized at Re = 500 
by perturbing simple Couette flow. We call the resulting 
unpatterned state uniform turbulence. Its spectrum is 
flat. Re is decreased quickly to 350 where a pattern forms 
with three distinct turbulent and laminar regions. The 
m = 3 spectral peak emerges. The selected wavelength of 
40 agrees closely with experiment 0, Q, . Re is kept 
at 350 long enough to show that this pattern is stable. 
The final flow at Re = 350 is visualized in Fig. ^ The 
pattern remains qualitatively the same through Re — 
320. E lms is systematically greater to the left of the band 
center. (Note that, due to the imposed tilt, there is no 
reflection symmetry in z' .) At Re = 310 the pattern loses 
one turbulent region, accompanied by the emergence of 
the m = 2 spectral peak. At Re = 300, a single turbulent 
region remains, and finally, at Re — 290, the flow reverts 
to simple Couette flow. 

Figure |3(b) shows states obtained by increasing Re 
starting from Re = 350. The steady three-banded pat- 
tern persists up through Re = 390. At Re = 400 and 
410 the pattern is no longer steady: bands are less well 
defined and laminar regions appear and disappear (see 
below). Uniform turbulence is obtained at Re = 420. 

We now present evidence that the patterns in Fig. 
represent three qualitatively different states. The banded 
state at Re = 350 is fundamentally spatially periodic. To 
support this we show in Fig.0]a simulation at Re = 350 in 
a domain whose length L z > is slowly increased. The pat- 
tern adjusts to keep the wavelength in the approximate 
range 35 — 65 by splitting the turbulent bands when they 
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FIG. 3: Space-time evolution of turbulent-laminar patterns 
in the domain L x i x L z > = 10 x 120, 8 = 24°. Time evolves 
upward with changes in Re indicated on the right. Grayscale 
or color plots: kinetic energy fluctuations E lms on a space- 
time grid. The same scale is used for all space-time plots, 
with E lma = in white. Insets: kinetic energy plotted over 
a time window T = 250 in a turbulent and laminar region. 
Left: Spectral peaks in the averaged spatial Fourier transform 
of kinetic energy with m — 3 (solid) and m = 2 (dotted), 
(long-dashed), and 1 (short-dashed), (a) States seen upon 
decreasing Re, from uniform turbulence at Re = 500, through 
various patterned states, ending in simple Couette flow at 
Re — 290. (b) States seen upon increasing Re, from the three- 
banded laminar-turbulent pattern at Re = 350 to uniform 
turbulence at Re = 420. 



grow too large. The instantaneous integrated energy pro- 
file E = J dx' dyE(x', y, z', t) is plotted at the final time. 
Between the turbulent bands, E does not reach zero and 
the flow, while basically laminar, differs from the simple 
Couette solution ye x . 

In sharp contrast, the single turbulent patch seen in 
Fig. |3|a) prior to return to laminar Couette flow is a 
localized state. Figure 0] shows that in a domain of in- 
creasing size at Re = 300 a single turbulent region of 
approximately fixed extent persists, independent of L z i. 
Moreover, E decays to zero exponentially as the flow 
approaches the simple Couette solution away from the 
patch. The localized states in our computations neces- 
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FIG. 4: Simulations at Re = 350, Re = 300, and Re = 410 
illustrating three distinct states: periodic, localized, intermit- 
tent. Space-time representation of E TnlB is as in Fig. |2l For 
Re — 350 and Re = 300 the domain length is increased from 
L z t = 50 to L z i — 140 in increments of 5. The integrated 
energy profile E(z') is shown at the final time. For Re = 410 
a single long simulation is shown for L z i = 40, accompanied 
by m = 1 (solid) and m — (dashed) spectral peaks. 

sarily take the form of bands when visualized in the x — z 
plane [e.g., Fig. Eld) below]. Isolated bands and spots are 
reported experimentally 0, Q, 0] near these values of 
Re, 

The third behavior is displayed by the intermittent 
state in Fig. [3] near the transition to uniform turbulence. 
Figure 0] shows a long simulation at Re = 410 in a do- 
main L z i = 40. The flow never stabilizes but instead 
quasi-laminar regions nucleate and disappear continually. 
The range of E rms in the space-time plot is noticeably 
smaller than for the stable patterns. Simulations with 
L z i = 60 show similar behavior. These states have been 
interpreted in [l3l [Til Ifo| as resulting from noise-driven 
competition between banded patterns at equal and op- 
posite angles. However, the intermittency is captured in 
our simulations, even though the competition between 
states of opposite angles is absent. 

We have increased and decreased Re gradually at 
L z i = 40 and L z > — 60 and find no hysteresis in any 
of the transitions between the turbulent states. 

We have explored regions of existence for various states 
as a function of Re, wavelength, and tilt. By varying L z : 
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FIG. 5: Turbulent-laminar patterns at minimum (6 — 15°) 
and maximum (8 = 66° ) angle for which they have been com- 
puted at Re = 350. Conventions as in Fig. Q 



ing two large lateral directions. Future studies of these 
states may shed light on the mechanisms responsible for 
laminar-turbulent patterns and for turbulent transition. 
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grant JR00WASTEQ). 



at 6 = 24°, Re — 350, we have determined that the 
minimum wavelength is 35 and the maximum is 65. For 
L z ' 30, uniform turbulence is obtained. For L z > > 70 
two bands of wavelength L z >/2 form (as in Fig.0J. This 
range of allowed wavelengths is nearly independent of Re 
wherever we have been able to compute banded states. 
Figure |3] shows a banded state at 9 — 15° and a lo- 
calized state at 9 = 66°, the minimum and maximum 
angles for which we have thus far obtained patterns for 
L 7 j = 120, Re — 350. These extreme states may not 
be stable without the imposed periodicity of the compu- 
tations. The sequence of states seen for increasing 9 at 
Re = 350 is qualitatively the same as that for decreasing 
Re at 9 = 24°. At 9 = 0° and 9 = 90° we do not find 
patterns, but only either uniform turbulence or simple 
Couette flow, with transition boundaries Re ~ 300 for 
9 = 0° and Re « 390 for 9 = 90°. Full details will be 
reported elsewhere. 

In past years minimal-flow-unit simulations have been 
used to great effect in understanding shear turbulence. 
We have shown that the same philosophy can be em- 
ployed in the study of large-scale structures formed in 
turbulent flows. Specifically, we have reported the first 
simulations of turbulent-laminar patterns in plane Cou- 
ette flow by numerically solving the Navier-Stokes equa- 
tions in domains with a single long direction. The other 
dimensions are just large enough to resolve the inter-plate 
distance and to contain an integer number of longitudinal 
vortex pairs or streaks. Thus we have demonstrated that 
the patterns are quasi one-dimensional and we have iden- 
tified what we believe to be near-minimal conditions nec- 
essary for their formation. Key is that the computational 
domain be tilted obliquely to the streamwise direction of 
the flow, otherwise no patterns are observed. We have 
found periodic, localized, and intermittent states where 
similar states are observed experimentally. We have ex- 
plored the patterns' dependence on Reynolds number and 
on imposed wavelength and tilt. The existence of lo- 
calized states in our simulations is particularly interest- 
ing because this suggests that the basic physics of iso- 
lated turbulent spots can be captured without simulat- 
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